Integrative multi-omics analysis of genomic, epigenomic, and metabolomics data leads to new insights for Attention-Deficit/Hyperactivity Disorder.

The evolving field of multi-omics combines data and provides methods for simultaneous analysis across several omics levels. Here, we integrated genomics (transmitted and non-transmitted polygenic scores [PGSs]), epigenomics, and metabolomics data in a multi-omics framework to identify biomarkers for Attention-Deficit/Hyperactivity Disorder (ADHD) and investigated the connections among the three omics levels. We first trained single- and next multi-omics models to differentiate between cases and controls in 596 twins (cases = 14.8%) from the Netherlands Twin Register (NTR) demonstrating reasonable in-sample prediction through cross-validation. The multi-omics model selected 30 PGSs, 143 CpGs, and 90 metabolites. We confirmed previous associations of ADHD with glucocorticoid exposure and the transmembrane protein family TMEM, show that the DNA methylation of the MAD1L1 gene associated with ADHD has a relation with parental smoking behavior, and present novel findings including associations between indirect genetic effects and CpGs of the STAP2 gene. However, out-of-sample prediction in NTR participants (N = 258, cases = 14.3%) and in a clinical sample (N = 145, cases = 51%) did not perform well (range misclassification was [0.40, 0.57]). The results highlighted connections between omics levels, with the strongest connections between non-transmitted PGSs, CpGs, and amino acid levels and show that multi-omics designs considering interrelated omics levels can help unravel the complex biology underlying ADHD.

Handedness and 23 Early Life Characteristics in 37,495 Dutch Twins.

In studies of singletons, a range of early-life characteristics have been reported to be associated with handedness, but some of these associations have failed to replicate. We examined associations between 23 early life characteristics with handedness in a large sample of 37,495 5-year-old twins. We considered three definitions of handedness: left-handedness (LH), mixed-handedness (MH), and non-right-handedness (NRH). Our main aim was to test whether the associations with sex, birth weight, gestational age, and season of birth - as reported in singletons - replicate in twins, and to examine twin-specific variables, including zygosity, chorionicity, birth order, and intertwin delivery time. Compared to previously published data from adults born as singletons (7.23%), the prevalence of NRH was higher in both twins (16.19%) and their parents (15.09%). In the twins, LH and NRH were associated with parents' LH. Male sex and lower gestational age were associated with NRH, and LH was associated with not being breastfed. MH was related to neurodevelopmental delays and higher externalizing problems later in childhood. Other previously reported associations were not replicated, and no twin-specific characteristics were related to handedness. These results emphasize the importance of considering multiple definitions of handedness and indicate a small number of replicated associations across studies.

Genetic and Environmental Effects on the Early Motor Development as a Function of Parental Educational Attainment.

INTRODUCTION: The contribution of genetic and environmental factors to individual differences in early motor development is still largely uncharted. This large-scale twin study establishes the genetic and environmental influences on the timing of motor milestones achievement, and it further tests whether the influences are moderated by parental education. METHODS: The twins came from families registered in the Netherlands Twin Register (NTR) from 1986 to 2016. In 30,256 complete twin pairs, mother-reported ages at which each twin was able to first-time roll from back to belly, sit unassisted, hands-and-knees crawl, stand up unaided, and walk independently were used to extract an early motor development factor. Parental education was dichotomized ("both parents with low/average education" vs "at least one parent with high education" with university degree as a threshold). RESULTS: Additive genetics explained 52% of the variance in motor development, the remaining 39% and 9% were explained by shared and nonshared environment separately. Mean age of achieving motor milestones tended to be higher in infants with high educated parents, and a moderation of parental education on the genetic and environmental variance in motor development was seen in female twins with larger heritability in the high educated parents group (64% vs 43%) paired to a lower shared environmental influence (28% vs 48%). Only 7%-8% of the variance was accounted for nonshared environmental factors, including measurement error. The pattern of results did not change when the degree of urbanicity, a correlate of parental education, was additionally considered. CONCLUSIONS: Genetic factors explain most of the individual differences in the timing of motor milestone achievement, but factors related to the shared home environment also play an important role in early motor development.

Integrative Multi-omics Analysis of Childhood Aggressive Behavior.

This study introduces and illustrates the potential of an integrated multi-omics approach in investigating the underlying biology of complex traits such as childhood aggressive behavior. In 645 twins (cases = 42%), we trained single- and integrative multi-omics models to identify biomarkers for subclinical aggression and investigated the connections among these biomarkers. Our data comprised transmitted and two non-transmitted polygenic scores (PGSs) for 15 traits, 78,772 CpGs, and 90 metabolites. The single-omics models selected 31 PGSs, 1614 CpGs, and 90 metabolites, and the multi-omics model comprised 44 PGSs, 746 CpGs, and 90 metabolites. The predictive accuracy for these models in the test (N = 277, cases = 42%) and independent clinical data (N = 142, cases = 45%) ranged from 43 to 57%. We observed strong connections between DNA methylation, amino acids, and parental non-transmitted PGSs for ADHD, Autism Spectrum Disorder, intelligence, smoking initiation, and self-reported health. Aggression-related omics traits link to known and novel risk factors, including inflammation, carcinogens, and smoking.

Changing genetic architecture of body mass index from infancy to early adulthood: an individual based pooled analysis of 25 twin cohorts.

BACKGROUND: Body mass index (BMI) shows strong continuity over childhood and adolescence and high childhood BMI is the strongest predictor of adult obesity. Genetic factors strongly contribute to this continuity, but it is still poorly known how their contribution changes over childhood and adolescence. Thus, we used the genetic twin design to estimate the genetic correlations of BMI from infancy to adulthood and compared them to the genetic correlations of height. METHODS: We pooled individual level data from 25 longitudinal twin cohorts including 38,530 complete twin pairs and having 283,766 longitudinal height and weight measures. The data were analyzed using Cholesky decomposition offering genetic and environmental correlations of BMI and height between all age combinations from 1 to 19 years of age. RESULTS: The genetic correlations of BMI and height were stronger than the trait correlations. For BMI, we found that genetic correlations decreased as the age between the assessments increased, a trend that was especially visible from early to middle childhood. In contrast, for height, the genetic correlations were strong between all ages. Age-to-age correlations between environmental factors shared by co-twins were found for BMI in early childhood but disappeared altogether by middle childhood. For height, shared environmental correlations persisted from infancy to adulthood. CONCLUSIONS: Our results suggest that the genes affecting BMI change over childhood and adolescence leading to decreasing age-to-age genetic correlations. This change is especially visible from early to middle childhood indicating that new genetic factors start to affect BMI in middle childhood. Identifying mediating pathways of these genetic factors can open possibilities for interventions, especially for those children with high genetic predisposition to adult obesity.

A genome-wide association study of total child psychiatric problems scores.

Substantial genetic correlations have been reported across psychiatric disorders and numerous cross-disorder genetic variants have been detected. To identify the genetic variants underlying general psychopathology in childhood, we performed a genome-wide association study using a total psychiatric problem score. We analyzed 6,844,199 common SNPs in 38,418 school-aged children from 20 population-based cohorts participating in the EAGLE consortium. The SNP heritability of total psychiatric problems was 5.4% (SE = 0.01) and two loci reached genome-wide significance: rs10767094 and rs202005905. We also observed an association of SBF2, a gene associated with neuroticism in previous GWAS, with total psychiatric problems. The genetic effects underlying the total score were shared with common psychiatric disorders only (attention-deficit/hyperactivity disorder, anxiety, depression, insomnia) (rG > 0.49), but not with autism or the less common adult disorders (schizophrenia, bipolar disorder, or eating disorders) (rG < 0.01). Importantly, the total psychiatric problem score also showed at least a moderate genetic correlation with intelligence, educational attainment, wellbeing, smoking, and body fat (rG > 0.29). The results suggest that many common genetic variants are associated with childhood psychiatric symptoms and related phenotypes in general instead of with specific symptoms. Further research is needed to establish causality and pleiotropic mechanisms between related traits.

DNA methylation in peripheral tissues and left-handedness.

Handedness has low heritability and epigenetic mechanisms have been proposed as an etiological mechanism. To examine this hypothesis, we performed an epigenome-wide association study of left-handedness. In a meta-analysis of 3914 adults of whole-blood DNA methylation, we observed that CpG sites located in proximity of handedness-associated genetic variants were more strongly associated with left-handedness than other CpG sites (P = 0.04), but did not identify any differentially methylated positions. In longitudinal analyses of DNA methylation in peripheral blood and buccal cells from children (N = 1737), we observed moderately stable associations across age (correlation range [0.355-0.578]), but inconsistent across tissues (correlation range [- 0.384 to 0.318]). We conclude that DNA methylation in peripheral tissues captures little of the variance in handedness. Future investigations should consider other more targeted sources of tissue, such as the brain.

Genetic factors explain a significant part of associations between adolescent well-being and the social environment.

Socio-environmental factors play an important role in adolescent well-being, but potential genetic contributions to these associations are rarely assessed. To address this gap in the literature, associations between well-being and family conflict and functioning, number of friends, friendship importance and satisfaction, and leisure time variables were studied in N = ~ 4700 twin pairs from the Netherlands Twin Register, us ing generalized estimating equations and twin-difference scores. When twin-difference scores indicated a role for genetic factors, we used bivariate genetic models to quantify genetic and environmental contributions to these associations. We identify significant associations between well-being and family functioning, family conflict, different leisure time activities, number of friends, and satisfaction with friendships. Additionally, we find evidence for large (73-91%) genetic influence on the associations between well-being and family conflict and functioning, leisure time sport/scouting clubs, and satisfaction with friendships. Finally, findings support the hypothesis of a causal association between well-being and family conflict and functioning. These findings have important implications for research into the social correlates of well-being in adolescence, as not taking genetic factors into account leads to overestimations of the influence of identified correlates and consequently to recommendations of these correlates as intervention targets.

Educational attainment of same-sex and opposite-sex dizygotic twins: An individual-level pooled study of 19 twin cohorts.

Comparing twins from same- and opposite-sex pairs can provide information on potential sex differences in a variety of outcomes, including socioeconomic-related outcomes such as educational attainment. It has been suggested that this design can be applied to examine the putative role of intrauterine exposure to testosterone for educational attainment, but the evidence is still disputed. Thus, we established an international database of twin data from 11 countries with 88,290 individual dizygotic twins born over 100 years and tested for differences between twins from same- and opposite-sex dizygotic pairs in educational attainment. Effect sizes with 95% confidence intervals (CI) were estimated by linear regression models after adjusting for birth year and twin study cohort. In contrast to the hypothesis, no difference was found in women (ß = -0.05 educational years, 95% CI -0.11, 0.02). However, men with a same-sex co-twin were slightly more educated than men having an opposite-sex co-twin (ß = 0.14 educational years, 95% CI 0.07, 0.21). No consistent differences in effect sizes were found between individual twin study cohorts representing Europe, the USA, and Australia or over the cohorts born during the 20th century, during which period the sex differences in education reversed favoring women in the latest birth cohorts. Further, no interaction was found with maternal or paternal education. Our results contradict the hypothesis that there would be differences in the intrauterine testosterone levels between same-sex and opposite-sex female twins affecting education. Our findings in men may point to social dynamics within same-sex twin pairs that may benefit men in their educational careers.

Identical twins carry a persistent epigenetic signature of early genome programming.

Monozygotic (MZ) twins and higher-order multiples arise when a zygote splits during pre-implantation stages of development. The mechanisms underpinning this event have remained a mystery. Because MZ twinning rarely runs in families, the leading hypothesis is that it occurs at random. Here, we show that MZ twinning is strongly associated with a stable DNA methylation signature in adult somatic tissues. This signature spans regions near telomeres and centromeres, Polycomb-repressed regions and heterochromatin, genes involved in cell-adhesion, WNT signaling, cell fate, and putative human metastable epialleles. Our study also demonstrates a never-anticipated corollary: because identical twins keep a lifelong molecular signature, we can retrospectively diagnose if a person was conceived as monozygotic twin.

Predicting Complex Traits and Exposures From Polygenic Scores and Blood and Buccal DNA Methylation Profiles.

We examined the performance of methylation scores (MS) and polygenic scores (PGS) for birth weight, BMI, prenatal maternal smoking exposure, and smoking status to assess the extent to which MS could predict these traits and exposures over and above the PGS in a multi-omics prediction model. MS may be seen as the epigenetic equivalent of PGS, but because of their dynamic nature and sensitivity of non-genetic exposures may add to complex trait prediction independently of PGS. MS and PGS were calculated based on genotype data and DNA-methylation data in blood samples from adults (Illumina 450 K; N = 2,431; mean age 35.6) and in buccal samples from children (Illumina EPIC; N = 1,128; mean age 9.6) from the Netherlands Twin Register. Weights to construct the scores were obtained from results of large epigenome-wide association studies (EWASs) based on whole blood or cord blood methylation data and genome-wide association studies (GWASs). In adults, MSs in blood predicted independently from PGSs, and outperformed PGSs for BMI, prenatal maternal smoking, and smoking status, but not for birth weight. The largest amount of variance explained by the multi-omics prediction model was for current vs. never smoking (54.6%) of which 54.4% was captured by the MS. The two predictors captured 16% of former vs. never smoking initiation variance (MS:15.5%, PGS: 0.5%), 17.7% of prenatal maternal smoking variance (MS:16.9%, PGS: 0.8%), 11.9% of BMI variance (MS: 6.4%, PGS 5.5%), and 1.9% of birth weight variance (MS: 0.4%, PGS: 1.5%). In children, MSs in buccal samples did not show independent predictive value. The largest amount of variance explained by the two predictors was for prenatal maternal smoking (2.6%), where the MSs contributed 1.5%. These results demonstrate that blood DNA MS in adults explain substantial variance in current smoking, large variance in former smoking, prenatal smoking, and BMI, but not in birth weight. Buccal cell DNA methylation scores have lower predictive value, which could be due to different tissues in the EWAS discovery studies and target sample, as well as to different ages. This study illustrates the value of combining polygenic scores with information from methylation data for complex traits and exposure prediction.

Large-scale collaboration in ENIGMA-EEG: A perspective on the meta-analytic approach to link neurological and psychiatric liability genes to electrophysiological brain activity.

BACKGROUND AND PURPOSE: The ENIGMA-EEG working group was established to enable large-scale international collaborations among cohorts that investigate the genetics of brain function measured with electroencephalography (EEG). In this perspective, we will discuss why analyzing the genetics of functional brain activity may be crucial for understanding how neurological and psychiatric liability genes affect the brain. METHODS: We summarize how we have performed our currently largest genome-wide association study of oscillatory brain activity in EEG recordings by meta-analyzing the results across five participating cohorts, resulting in the first genome-wide significant hits for oscillatory brain function located in/near genes that were previously associated with psychiatric disorders. We describe how we have tackled methodological issues surrounding genetic meta-analysis of EEG features. We discuss the importance of harmonizing EEG signal processing, cleaning, and feature extraction. Finally, we explain our selection of EEG features currently being investigated, including the temporal dynamics of oscillations and the connectivity network based on synchronization of oscillations. RESULTS: We present data that show how to perform systematic quality control and evaluate how choices in reference electrode and montage affect individual differences in EEG parameters. CONCLUSION: The long list of potential challenges to our large-scale meta-analytic approach requires extensive effort and organization between participating cohorts; however, our perspective shows that these challenges are surmountable. Our perspective argues that elucidating the genetic of EEG oscillatory activity is a worthwhile effort in order to elucidate the pathway from gene to disease liability.

Associations of sleep with psychological problems and well-being in adolescence: causality or common genetic predispositions?

BACKGROUND: Whereas short and problematic sleep are associated with psychological problems in adolescence, causality remains to be elucidated. This study therefore utilized the discordant monozygotic cotwin design and cross-lagged models to investigate how short and problematic sleep affect psychological functioning. METHODS: Adolescent twins (N = 12,803, 13-20 years, 42% male) completed questionnaires on sleep and psychological functioning repeatedly over a two-year interval. Monozygotic twin pairs were classified as concordant or discordant for sleep duration and trouble sleeping. Resulting subgroups were compared regarding internalizing problems, externalizing problems, and subjective well-being. RESULTS: Cross-sectional analyses indicated associations of worse psychological functioning with both short sleep and problematic sleep, and cross-lagged models indicate bidirectional associations. Longitudinal analyses showed that an increase in sleep problems experienced selectively by one individual of an identical twin pair was accompanied by an increase of 52% in internalizing problem scores and 25% in externalizing problem scores. These changes were significantly different from the within-subject changes in cotwins with unchanged sleep quality (respectively, 3% increase and 5% decrease). Psychological functioning did, however, not worsen with decreasing sleep duration. CONCLUSIONS: The findings suggest that sleep quality, rather than sleep duration, should be the primary target for prevention and intervention, with possible effect on psychological functioning in adolescents.

Early-life antibiotic use and risk of attention-deficit hyperactivity disorder and autism spectrum disorder: results of a discordant twin study.

BACKGROUND: Development of the gut-brain axis in early life may be disturbed by antibiotic use. It has been hypothesized that this disturbance may contribute to development of neurodevelopmental disorders, including autism spectrum disorder and attention-deficit hyperactivity disorder. We aimed to assess the association between antibiotic use in early life and the risk of developing attention-deficit hyperactivity disorder or autism spectrum disorder, while controlling for shared genetic and environmental factors in a discordant twin design. METHODS: We conducted a cohort study in twins (7-12 years; 25 781 twins) from the Netherlands Twin Register (NTR) and a replication study in the Childhood and Adolescent Twin Study in Sweden (CATSS; 7946 9-year-old twins). Antibiotic use was recorded before age 2 years. Attention-deficit hyperactivity disorder and autism spectrum disorder were parent-reported in the Netherlands Twin Register and register-based in the Childhood and Adolescent Twin Study in Sweden. RESULTS: Early-life antibiotic use was associated with increased risk of attention-deficit hyperactivity disorder development [pooled odds ratio (OR) 1.10, 95% confidence interval (CI) 1.02-1.17] and autism spectrum disorder (pooled OR 1.15, 95% CI 1.06-1.25) in a case-control design. When restricting to monozygotic twin pairs discordant for the outcome, associations disappeared for both disorders in both cohorts (attention-deficit hyperactivity disorder OR 0.90, 95% CI 0.48-1.69 and OR 0.80, 95% CI 0.37-1.76, and autism spectrum disorder OR 0.66, 95% CI 0.38-1.16 and OR 0.29, 95% CI 0.02-4.50, respectively). CONCLUSIONS: Our findings suggest that the association between early-life antibiotic use and risk of attention-deficit hyperactivity and autism spectrum disorder may be confounded by shared familial environment and genetics.

Urinary Amine and Organic Acid Metabolites Evaluated as Markers for Childhood Aggression: The ACTION Biomarker Study.

Biomarkers are of interest as potential diagnostic and predictive instruments in personalized medicine. We present the first urinary metabolomics biomarker study of childhood aggression. We aim to examine the association of urinary metabolites and neurotransmitter ratios involved in key metabolic and neurotransmitter pathways in a large cohort of twins (N = 1,347) and clinic-referred children (N = 183) with an average age of 9.7 years. This study is part of ACTION (Aggression in Children: Unraveling gene-environment interplay to inform Treatment and InterventiON strategies), in which we developed a standardized protocol for large-scale collection of urine samples in children. Our analytical design consisted of three phases: a discovery phase in twins scoring low or high on aggression (N = 783); a replication phase in twin pairs discordant for aggression (N = 378); and a validation phase in clinical cases and matched twin controls (N = 367). In the discovery phase, 6 biomarkers were significantly associated with childhood aggression, of which the association of O-phosphoserine (ß = 0.36; SE = 0.09; p = 0.004), and gamma-L-glutamyl-L-alanine (ß = 0.32; SE = 0.09; p = 0.01) remained significant after multiple testing. Although non-significant, the directions of effect were congruent between the discovery and replication analyses for six biomarkers and two neurotransmitter ratios and the concentrations of 6 amines differed between low and high aggressive twins. In the validation analyses, the top biomarkers and neurotransmitter ratios, with congruent directions of effect, showed no significant associations with childhood aggression. We find suggestive evidence for associations of childhood aggression with metabolic dysregulation of neurotransmission, oxidative stress, and energy metabolism. Although replication is required, our findings provide starting points to investigate causal and pleiotropic effects of these dysregulations on childhood aggression.

Content, diagnostic, correlational, and genetic similarities between common measures of childhood aggressive behaviors and related psychiatric traits.

BACKGROUND: Given the role of childhood aggressive behavior (AGG) in everyday child development, precise and accurate measurement is critical in clinical practice and research. This study aims to quantify agreement among widely used measures of childhood AGG regarding item content, clinical concordance, correlation, and underlying genetic construct. METHODS: We analyzed data from 1254 Dutch twin pairs (age 8-10 years, 51.1% boys) from a general population sample for whom both parents completed the A-TAC, CBCL, and SDQ at the same occasion. RESULTS: There was substantial variation in item content among AGG measures, ranging from .00 (i.e., mutually exclusive) to .50 (moderate agreement). Clinical concordance (i.e., do the same children score above a clinical threshold among AGG measures) was very weak to moderate with estimates ranging between .01 and .43 for mother-reports and between .12 and .42 for father-reports. Correlations among scales were weak to strong, ranging from .32 to .70 for mother-reports and from .32 to .64 for father-reports. We found weak to very strong genetic correlations among the measures, with estimates between .65 and .84 for mother-reports and between .30 and .87 for father-reports. CONCLUSIONS: Our results demonstrated that degree of agreement between measures of AGG depends on the type (i.e., item content, clinical concordance, correlation, genetic correlation) of agreement considered. Because agreement was higher for correlations compared to clinical concordance (i.e., above or below a clinical cutoff), we propose the use of continuous scores to assess AGG, especially for combining data with different measures. Although item content can be different and agreement among observed measures may not be high, the genetic correlations indicate that the underlying genetic liability for childhood AGG is consistent across measures.

Harmonizing behavioral outcomes across studies, raters, and countries: application to the genetic analysis of aggression in the ACTION Consortium.

BACKGROUND: Aggression in children has genetic and environmental causes. Studies of aggression can pool existing datasets to include more complex models of social effects. Such analyses require large datasets with harmonized outcome measures. Here, we made use of a reference panel for phenotype data to harmonize multiple aggression measures in school-aged children to jointly analyze data from five large twin cohorts. METHODS: Individual level aggression data on 86,559 children (42,468 twin pairs) were available in five European twin cohorts measured by different instruments. A phenotypic reference panel was collected which enabled a model-based phenotype harmonization approach. A bi-factor integration model in the integrative data analysis framework was developed to model aggression across studies while adjusting for rater, age, and sex. Finally, harmonized aggression scores were analyzed to estimate contributions of genes, environment, and social interaction to aggression. The large sample size allowed adequate power to test for sibling interaction effects, with unique dynamics permitted for opposite-sex twins. RESULTS: The best-fitting model found a high level of overall heritability of aggression (~60%). Different heritability rates of aggression across sex were marginally significant, with heritability estimates in boys of ~64% and ~58% in girls. Sibling interaction effects were only significant in the opposite-sex twin pairs: the interaction effect of males on their female co-twin differed from the effect of females on their male co-twin. An aggressive female had a positive effect on male co-twin aggression, whereas more aggression in males had a negative influence on a female co-twin. CONCLUSIONS: Opposite-sex twins displayed unique social dynamics of aggressive behaviors in a joint analysis of a large, multinational dataset. The integrative data analysis framework, applied in combination with a reference panel, has the potential to elucidate broad, generalizable results in the investigation of common psychological traits in children.

Out of Control: Examining the Association Between Family Conflict and Self-Control in Adolescence in a Genetically Sensitive Design.

OBJECTIVE: Family conflict is associated with low self-control in adolescence. Thus far research about the direction of this association is inconclusive. In this study, we sort out whether this association reflects a causal effect or whether it is explained by a common underlying cause, including genetic factors. METHOD: In twin data, we fitted a series of causal models, and compared models for the association of family conflict and self-control including reciprocal causation, unidirectional causation from family conflict to low self-control, unidirectional causation from low self-control to family conflict, and common genetic susceptibility. We included data from a large sample of twins aged 14 years (N = 9,173), all enrolled in the Netherlands Twin Register. RESULTS: The results suggested a unidirectional pathway model in which family conflict leads to low self-control in adolescence, with genetic factors also playing a role in explaining the association. CONCLUSION: Adolescents experiencing family conflict are at risk for showing hampered self-control capacities, with family conflict being a robust predictor of low self-control through common genetic factors but also through direct causal influences.

Parental Age and Offspring Childhood Mental Health: A Multi-Cohort, Population-Based Investigation.

To examine the contributions of maternal and paternal age on offspring externalizing and internalizing problems, this study analyzed problem behaviors at age 10-12 years from four Dutch population-based cohorts (N = 32,892) by a multiple informant design. Bayesian evidence synthesis was used to combine results across cohorts with 50% of the data analyzed for discovery and 50% for confirmation. There was evidence of a robust negative linear relation between parental age and externalizing problems as reported by parents. In teacher-reports, this relation was largely explained by parental socio-economic status. Parental age had limited to no association with internalizing problems. Thus, in this large population-based study, either a beneficial or no effect of advanced parenthood on child problem behavior was observed.

DNA Methylation Signatures of Breastfeeding in Buccal Cells Collected in Mid-Childhood.

Breastfeeding has long-term benefits for children that may be mediated via the epigenome. This pathway has been hypothesized, but the number of empirical studies in humans is small and mostly done by using peripheral blood as the DNA source. We performed an epigenome-wide association study (EWAS) in buccal cells collected around age nine (mean = 9.5) from 1006 twins recruited by the Netherlands Twin Register (NTR). An age-stratified analysis examined if effects attenuate with age (median split at 10 years; n<10 = 517, mean age = 7.9; n>10 = 489, mean age = 11.2). We performed replication analyses in two independent cohorts from the NTR (buccal cells) and the Avon Longitudinal Study of Parents and Children (ALSPAC) (peripheral blood), and we tested loci previously associated with breastfeeding in epigenetic studies. Genome-wide DNA methylation was assessed with the Illumina Infinium MethylationEPIC BeadChip (Illumina, San Diego, CA, USA) in the NTR and with the HumanMethylation450 Bead Chip in the ALSPAC. The duration of breastfeeding was dichotomized ('never' vs. 'ever'). In the total sample, no robustly associated epigenome-wide significant CpGs were identified (α = 6.34 × 10-8). In the sub-group of children younger than 10 years, four significant CpGs were associated with breastfeeding after adjusting for child and maternal characteristics. In children older than 10 years, methylation differences at these CpGs were smaller and non-significant. The findings did not replicate in the NTR sample (n = 98; mean age = 7.5 years), and no nearby sites were associated with breastfeeding in the ALSPAC study (n = 938; mean age = 7.4). Of the CpG sites previously reported in the literature, three were associated with breastfeeding in children younger than 10 years, thus showing that these CpGs are associated with breastfeeding in buccal and blood cells. Our study is the first to show that breastfeeding is associated with epigenetic variation in buccal cells in children. Further studies are needed to investigate if methylation differences at these loci are caused by breastfeeding or by other unmeasured confounders, as well as what mechanism drives changes in associations with age.